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Abstract. The graph isomorphism problem has a long history in math- 
ematics and computer science, with applications in computational chem- 
istry and biology, and it is believed to be neither solvable in polynomial 
time nor NP-complete. E. Luks proposed in 1982 the best algorithm so 
far for the solution of this problem, which moreover runs in polynomial 
time if an upper bound for the degrees of the nodes in the graphs is taken 
as a constant. Unfortunately, Luks' algorithm is purely theoretical, very 
difficult to use in practice, and, in particular, we have not been able to 
find any implementation of it in the literature. The main goal of this pa- 
per is to present an efficient implementation of this algorithm for ternary 
graphs in the SAGE system, as well as an adaptation to fully resolved 
rooted phylogenetic networks on a given set of taxa. 



1 Introduction 

The graph isomorphism problem — the problem of deciding whether two finite 
graphs are isomorphic or not — has a long history in mathematics and computer 
science. It is one of the most important decision problems for which the compu- 
tational complexity is not known yet |11|16 : , and it is believed to be neither in 
P nor NP-complete. Its uncertain status has led to the definition and study of 
the class GI of the decision problems that are polynomial-time Turing reducible 
to it [13]. 

The best current theoretical algorithm for the solution of the graph iso- 
morphism problem is due to Eugene Luks |18|19j . His algorithm runs in time 
20(yniog(n)\ where n is the number of nodes in the graphs, and it relies on the 
classification of finite simple groups [6 j. Without using this major result, whose 
proof consists of tens of thousands of pages in several hundred journal articles 
and for which a unified, simplified and revised version is still in progress, Babai 
and Luks [3] gave an algorithm that runs in 2°(^ log ( n ) ) time. 

The graph isomorphism problem has many practical applications outside 
mathematics. Our interest in it stems from its applications in computational 
biology and bioinformatics. Graphs are ubiquitous in biology as models of dif- 
ferent complex systems: molecular structures, phylogenetic trees and networks, 



metabolical pathways, protein-protein interaction (PPI) networks, gene expres- 
sion networks, etc. |2|14|15)20] . In all these fields, the comparison of graphs is an 
important computational problem. For instance, a protein's function is closely 
related to its three-dimensional shape, which at an intermediate level of detail 
is modeled by a contact graph, and thus the comparison of such contact graphs 
is a key tool in the prediction of proteins' function [12 23 ; the comparison of 
phylogenies is acknowledged to be one of the main problems in phylogenetics 
|7|15|22| : and the comparison of metabolic pathways and PPI networks is cur- 
rently one of the hot tools in the study of evolution |8|17j . The most basic graph 
comparison problem in all these contexts is, of course, the detection of equalities. 
Without being able to decide efficiently whether two such graphs are isomorphic 
or not, we cannot expect to define metrics, matchings or alignments that can be 
computed efficiently. 

So, it is important to design efficient algorithms to test if two graphs are 
isomorphic, at least for some special classes of graphs with practical applica- 
tions. The aforementioned algorithm by Luks [19] is fixed-parameter tractable, 
in the sense that it runs in polynomial time if we consider an upper bound d for 
the degrees of the graphs under comparison as a constant. Such an algorithm, 
even for ternary graphs (graphs with all their nodes of degree at most 3), has 
relevant applications in computational biology: for instance, in phylogenetics, 
in the comparison of arbitrary fully resolved rooted phylogenetic networks, and 
in the comparison of split networks [15]; in both types of networks, all nodes 
have total degrees at most 3. For these networks, no specific polynomial-time 
isomorphism test has been devised yet, and the authors working with them sim- 
ply quote Luks' result when they need to state that their isomorphism can be 
decided in polynomial time; see, for instance, [15] p. 168] or [4j §V]. 

The problem with Luks' algorithm is that it is purely theoretical, very dif- 
ficult to use in practice, and, in particular, we have not been able to find any 
implementation of it in the literature. The main goal of this paper is, then, to 
present an efficient implementation of this algorithm in the SAGE system. We 
also present an implementation of an adaptation of this algorithm that solves 
the isomorphism problem for fully resolved rooted phylogenetic networks. This 
paper is based on the first author's Master Thesis pQ, where we refer the reader 
for more details. 

2 Luks' algorithm for ternary graphs 

In this section we explain Luks' algorithm for ternary graphs, that is, for con- 
nected graphs with all their nodes of degree at most 3. We shall omit all proofs, 
which can be found in the original paper by Luks [19] or, with more detail, in 
the first author's Master Thesis pQ. 

We begin by reducing the graph isomorphism problem to a problem about 
isomorphisms of groups. Given two ternary graphs X\ = (Vi,Ei) and X 2 = 
(1/2,^2)? an d two edges e\ G E\ and e2 G £"2, consider the ternary graph X = 
BuildX(Xi, X 2 , ei, e 2 ) obtained by splitting e\ and e 2 by means of new nodes V\ 



and V2 , respectively, and then connecting these new nodes v\ and V2 by means of 
a new edge e. Then, X\ and X2 are isomorphic if, and only if, given any e\ G Ei, 
there exist some 62 G E2 and some automorphism a of X such that cr(i?i) = V2, 
and in particular such that a(e) = e. And then, if any such automorphism of X 
does exist, then any set of generators of Aut e (X), the group of automorphisms 
of G that fixes e, will contain some. This provides the following algorithm. 



Algorithm 1: Isomorphism of ternary graphs 



Data: X\ , X2 ternary graphs 
Result: Test if X\ and X2 are isomorphic 

1 begin 

2 ei G E{X 1 ) 

3 for e 2 G £(X 2 ) do 

4 X <- BuildX(Xi,X 2 ,ei,e 2 ) 

5 5 <(— a set of generators of Aut e (X) 

6 for a G S do 

7 if cr(vi) v 2 then 

8 return True 

9 return False 



So, we are reduced to compute a set of generators of Aut e (X), for a ternary 
graph X and any e G E(X). This group is determined through a natural sequence 
of successive "approximations" Aut e (X r ), where each X r is the subgraph con- 
sisting of all nodes and all edges of X which appear in paths of length < r 
containing e. More formally, if e = (a, 6), let 

X 1 = ({a,b},{(a,b)}) 

and, for every r > 2, 

V(X r ) = {ye V(X) I 3x G V(X r _ 1 ) such that G E'(X)} 

£;(X r ) = G I 3x G V(X r _i) such that (x,y) G ^(AT)} 

There exist natural homomorphisms 

7r r : Aut e (X r+ i) — >> Aut e (X r ) 

defined by the restriction of the automorphisms. These homomorphisms allow us 
to construct recursively a generating set for Aut e (X r+ i) from one for Aut e (X r ) 
by solving, for every r > 1, the following two problems: 

(I) Compute a set /C r of generators of ker7r r . 
(II) Compute a set <S r+ i of generators of 7r r (Aut e (X r+ i)). 



Indeed, if S' r+1 is any set of pre-images of <S r+ i in Aut e (X r+ i), then JC r U S' r+1 
generates Aut e (X r+ i). 

So, if we know how to solve Problems I and II, the algorithm to compute a 
set of generators of Aut e (X) is the following. In it, and in the sequel, given two 
elements a, b of a set F, we denote by (a b) G Sym(F) the transposition a o b. 



Algorithm 2: A set of generators of Aut e (X) 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 



Data: A ternary graph X and an edge e 
Result: A set of generators of Aut e (X) 
begin 

Compute the sequence of subgraphs X\ C • • • C X]\ 
Aut e = {(ab)} 
while r < N do 

K <(— a set of generators of ker 7r r 
S ^— a set of generators of 7r r (Aut e (Xr+i)) 
S' <(— a set of pre-images of S in Aut e (X r+ i) 
Aut e = S' U K 
r = r + 1 

return Aut e 



X as before 



Set now V r+1 = V(X r+1 ) \ V(X r ) and A r = {a C V r \ < \a\ < 3}, and 
consider the mapping 

/ : V r+1 A r 

defined by f(v) = {w G V(X r ) \ (v,w) G E(X)}; we call f(v) the neighbor set 
of v. The following result solves Problem I. 

Proposition 1. JC r = {(u v) G Sym(V r+ i) | ^ v, f(u) = /(v)}. 

Moreover, the fact that each ker7r r is generated by transpositions implies, by 
induction, the following result, which will be useful later. 

Proposition 2. For each r, Aut e (X r ) is a 2-group. 

So, our graph isomorphism problem is reduced to solve Problem II. Now, 
consider the following three subsets of A r : 

A! r = {{v u v 2 } G A r | (v u v 2 ) G E(X r+l )} 

A r ,i = {a G A r | a = f(v) for some unique v G K+i} 

A r ,2 = {a G A r | a = = /(^2) for some v\ ^ v 2 } 

We have the following result. 



Proposition 3. 7r r (Aut e (X r+ i)) is the set of those a G Aut e (X r ) that stabilize 
the sets A rj ±, A r $, and A' r . 

Now, set B r = V(X r -i) U A r and extend the action of Aut e (X r ) on V(X r ) 
to B r in the natural way: if a G A r , then a (a) = {o~(w) \ w G a}. Color each 
element of B r with one of five colors that distinguish, on the one hand, whether 
or not it is in A f r , and, on the other hand, whether it is in A ri i, or A r ^-, or 
neither. Only five colors are needed, since A' r n A r ^ = 0. 

By the previous proposition, a G 7r r (Aut e (X r+ i)) if and only if a preserves 
these colors in A r . Therefore, the isomorphism problem for ternary graphs is 
polynomial-time reducible to the following color automorphism problem (taking 
G = Aut e (X r ), B = A = A r , a = Id). 

Problem 1. Given 

— A set of generators for a 2-subgroup G of the symmetric group Sym(A) of a 
colored set A 

— A G-stable subset B C A 

— A permutation a G Sym(A) 

find Cb{o~G) 1 where, for every T C Sym(A) 

Cb(T) = {r G T | r preserves the color of every b G B} 

Now, the following three lemmas are the basis of Algorithm 3 that solves 
Problem 1 in polynomial-time, and thus it completes the graph isomorphism 
test we were looking for. 

Lemma 1. Let G be a subgroup of Syn^A), a G Syu^A) and B a G -stable 
subset of A such that Cb(o~G) is not empty. Then, Cb(G) is a subgroup of G 
and Cb(o~G) it is a left coset of the subgroup Cb(G). 

Lemma 2 (Furst-Hopcroft-Luks [9]). Given a set of generators of a sub- 
group G of a symmetric group, one can compute in polynomial time a set of 
generators of any subgroup of G that is known to have polynomially bounded 
index in G and for which a polynomial-time membership test is available. 

Lemma 3 (Luks [19j). Given a set of generators of a subgroup G of a sym- 
metric group and a G-orbit B, one can determine in polynomial time a minimal 
G -block system in B. 

Theorem 1. Algorithm 1 solves the ternary graph isomorphism in time 0(n 10 ), 
with n the size of the input graphs. 

Proof. The cost of computing BuildX applied to a pair of graphs of size n and 
a pair of edges is in 0(n) time. For every r, a set of generators of ker7r r can 
be computed as explained in Lemma 1 in 0(n 2 ) time, a set of generators of 
7r r (Aut e (X r+ i)) can be computed, by reducing this computation to Problem 1 



Algorithm 3: C B (crG) 

Data: A coset aG C Sym(A), where A is a colored set and G is a 2-group, and 

a G- stable subset B of A. 
Result: C B (aG) 



l begin 


2 


case D — {h} 


3 




if a{b) ~ b then 


4 




| Cb(o-G) = aG 


5 




else 


6 




L c B = ® 


7 


case G is intransitive on B 


8 




Let a nontrivial orbit 


9 




B 2 = B \ B x 


10 




_ Cb(o-G) = Cb 2 (C Bi (vG)) 


11 


case G is transitive on B 


12 




Let {Bi, B2} a minimal G-block system 


13 




Find the subgroup H of G that stabilizes B\ 


14 




Let r e G\H 


15 




_ Cb(o-G) = Cb 2 (Cb 1 (vH))UCb 2 (C Bi (vtH)) 


16 


return Cb(o~G) 



and then using Algorithm 3, in 0(n 8 ) time, and a set of preimages in Aut e (X r+ i) 
of the latter can be computed in 0(n 2 ) time. Since there are at most 0(n) indices 
r, we conclude that the cost of computing a set of generators of Aut e (X) by 
means of Algorithm 2 is in 0(n 9 ). Since Algorithm 1 calls Algorithm 2 0{n) 
times, the final cost of Algorithm 1 is in 0(n 10 ). 

3 Implementation 
3.1 Improvements 

In our implementation we have improved the efficiency of Luks' algorithm by: 

(a) Reducing the size of the set A r . 

(b) Representing the groups by means of smooth generating sequences. 

(c) Precomputing the blocks. 

(d) Running initial tests, to avoid trivial cases. 

(e) Removing, for every r, the permutations that do not swap the two "parts" 
Xl and X 2 of X. 

With these improvements, the cost of the computation decreases to 0(n 4 ). We 
explain now in some detail the improvements (a)-(c), which were inspired by 
[TO] . For more details, see pQ. 

As far as improvement (a) goes, we have been able to remove the triplets from 
B ri by replacing each node v with a 3-elements neighbor set by a triangle with 



Fig. 1. Replacing the triplets in the neighbor sets 



nodes at "level" r + 1 and labeled edges: cf. Figure 3.1. In this way, we can replace 
the graph X by a new graph X with some edges labeled. The automorphisms of 
X must preserve labels, and therefore the computation of Aut e (X) is the same 
as that of Aut e (X), except for the following facts: B r needs only to include the 
subsets of V r of size 1 or 2; we split A' r into 

A r,a = {{^1^2} G A r I (vi,v 2 ) e E(X r+1 ) is unlabeled}, 
A r,b = {{^1^2} e A r I (v!,v 2 ) e E(X r+1 ) is labeled}; 

and we modify the set of colors on B r to distinguish, on the one hand, whether 
or not an element is in A' r a , or A' r 6 , or neither, and, on the other hand, whether 
it is in A r i, or A r ^, or neither. 

As to (b), we represent 2-groups in a way that makes easier several key 
computations. 

Definition 1. Let G be a 2-group generated by {#i, . . . ,#&}. The sequence (#1, . . . , 
is a smooth generating sequence (SGS, for short) for G if [{git • . . ,gi) : (^1, . . . ,gi-i}] 
2, for i = 1, . . . , k. 

SGS are preserved by homomorphisms and liftings, and if we know a SGS 
for a 2-group G, then it is easy construct an SGS for a subgroup H of index 2 

m- 

Lemma 4. Let G be a 2-group, (#i, . . . , gk) a SGS for it, and H a subgroup of 
index 2. Let j = mm{i | gi £ H} and set, for i = 1, . . . , k, 

o = / 9i if 9i^H 

Then . . . , fa) is a SGS for H, and this sequence is computed in O(k) times 
the required time of a membership test for H. 

Let us consider finally improvement (c). As we have seen in the proof of 
Theorem 1, the most expensive part of Luks' algorithm is the recursive calls 
performed by Algorithm 3. The task carried out by this algorithm can be reor- 
ganized so as to limit the number of different blocks visited. These blocks form 
a tree that is precomputed and guides the recursion. 



Definition 2. Let G be a 2-group acting on a stable subset B of a colored set A. 
We call a binary tree T a structure tree for B with respect to G, T = T(£?, G), 
if the set of leaves of T is B, and the action of any a G G on B can be lifted to 
an automorphism of T . 

Let Q be a fixed color of A (in our application, it will be the color of the 
elements of A r that do not belong to A! r U A r ^ U A r>2 ). A node B of T = T(B, G) 
will be called active if B D Q c ^ 0. An active node 5 is facile if G is intransitive 
on B and the latter has exactly one active child. Let A (B) denote the closest 
non facile descendant of B in T. 

We can precompute the entire structure tree T(B, G) for a given pair (£?, G) 
using Algorithm 4 below. With this algorithm, we can construct every structure 
tree T(B ri G r ) in 0{n 2 ) time, and we can compute the mapping A and the set 
of active nodes in 0(n log n) time. Using structure trees to guide the recursion, 
the cost of computing Cb(o~G) decreases to 0(n 4 ). 



Algorithm 4: T(B, G) 



Data: B, G 

Result: T — T(B, G) 

1 begin 

2 Let the root of T be B 

3 if \B\ 1 then 

4 return 

5 Find the orbits of G in B 

6 if G is transitive then 

7 Find a minimal block system {Bl, Br} for G on B 

8 Find the subgroup H of G that stabilizes Bz, 

9 Find r e G\H 

10 return T = Tree(B L , H) U r(Tree(B L , H) (linked to the new root B ) 

n else 

12 Partition B into two nontrivial G-stable subsets Bl, Br 

13 return T = T(B L , G) U T(B R , G) (linked to the new root B) 



3.2 Implementation details 

We have implemented Luks' algorithm using the language Python and some 
specific SAGE libraries for handling graphs and groups of permutations. 

Besides the obvious classes and functions necessary to implement the algo- 
rithm, we use a new class, called Node, for the structure tree. This class has four 
attributes: 

— Node, the content of the node. 



— Left, the left child. 

— Right, the right child. 

— Parent, the parent. 

When we build a new Node without some attributes, they will be empty arrays. 
This class has the necessary functions to modify its attributes, as well as the 
following functions: 

— IsLeaf (), to know whether the node is a leaf. 

— IstransitiveO , to know whether G is transitive on the node. 

— Isactive(Q), to know whether the node is active w.r.t. Q. 

— Isfacile(Q), to know whether the node is facile w.r.t. Q. 

— Delta(Q), to know the nearest non facile descendant of the node. 

To test whether two graphs are isomorphic or not, we can use two different 
functions, Isomorphism and Isomorphism2. Both functions answer the question 
whether the input graphs are isomorphic, but moreover the first function returns 
the whole group of automorphisms that fix the distinguished edge e of the graph 
X = BuildXpTi, X 2: ei, e 2 ) : while the second one only returns the subgroup of 
those automorphisms of X that swap the parts X\ and X 2 . 

Finally, we have adapted Luks' algorithm to test the isomorphism of fully re- 
solved rooted phylogenetic networks on a given set of taxa. A rooted phylogenetic 
network, on a set of taxa S is a rooted, directed, acyclic graph with its leaves 
bijectively labeled in S. These graphs are used as explicit models of evolution- 
ary histories that, besides mutations, include reticulate evolutionary events like 
genetic recombinations, lateral gene transfers or hybridizations. An evolutionary 
network is fully resolved , or binary, when, for every node v in it, the ordered 
pair (di n (v),d ou t(v)) is either (0,2) (the root), (1,0) (the leaves), (1,2) (the 
tree nodes) or (2, 1) (the reticulate nodes). Two evolutionary networks on S are 
isomorphic when they are isomorphic as directed graphs and the isomorphism 
preserves the leaves' labels. For more on phylogenetic networks, see [15]. 

The main difference between the case of rooted phylogenetic networks and 
the general case is that, in the former, the isomorphisms map the root to the 
root, and therefore the graph X = BuildX(Xi, X 2 ) can be simply obtained by 
connecting the roots by an edge e. In particular, Algorithm 1 needs not to call 
0(n) times Algorithm 2, but only once, and the resulting cost is then in 0(n 3 ). 

To test the isomorphism of fully resolved rooted phylogenetic networks we 
have defined the function IsomorphismPhilo, which has the following parame- 
ters: 

— XI, the first phylogenetic network. 

— X2, the second phylogenetic network. 

— n, the number of nodes of XI. 

— dicl, the dictionary with the labels of the nodes of the first graph. 

— dic2, the dictionary with the labels of the nodes of the second graph. 

— rl, the root of the first graph. 

— r2, the root of the second graph. 



IsomorphismPhilo accepts networks with inner taxa (with internal labeled 
nodes) as well as multilabeled networks (where different nodes can have the same 
label), although these features are not used (yet) in the phylogenetic networks 
literature. For the moment, the algorithm only accepts phylogenetic networks 
created with SAGE's function Graph, but in the near future we plan to adapt it 
so that it accepts networks described in the format eNewick [5]. 

We have performed several tests on our implementation, which we report in 
the next subsection. In the first and the third tests, we have used SAGE functions 
to generate random graphs and we have then made these graphs connected and 
trivalent, by randomly adding or removing edges when it was necessary. In the 
second test, we have used the function graphs .DegreeSequence that, given a 
sequence of degrees, returns a graph whose nodes have this sequence of degrees, 
if some exists. This has allowed us to know the probability that the graphs under 
comparison were isomorphic. 

To perform the tests on phylogenetic networks we defined some functions to 
create random fully resolved rooted phylogenetic networks: 

— createDic (nodes, n), which, given a set of nodes, returns the following dic- 



— createDic2 (nodes ,n) , as above, but the selected nodes are labeled "a" or 
"b" equiprobably. 

— leaves (XI), returns the set of nodes of the network XI. 

Then, we have developed an algorithm RandomTree(n) that returns a fully re- 
solved rooted phylogenetic network with n nodes. By default, each internal node 
has a probability of 0.5 of being hybrid, but it can be changed by simply changing 
the probability parameter from 0.5 to the desired probability. 

The documentation of the whole module can be found in http : //www . alumnos . 
unican . es/ aam35/ sage-epydoc/ index . html and the code in http : //www . alumnos . 
unican . es/aam35/IsoTriGraph . py 



3.3 Tests 

The first two examples show that the code works correctly, and then the tests 
prove that it runs in a reasonable time, comparable to the speed of the own 
SAGE algorithm to test isomorphisms of graphs. 

Example 1. Consider the graphs X\ and X2 depicted in Fig. 2. These graphs 
are created in SAGE with: 



sage: X3=Graph( [(1 , 7), (1, 10), (2, 3), (2, 4), (3, 4), (4, 9), 

(5, 6), (6, 8), (7, 8), (7, 9), (8, 9)]) 
sage: X4=Graph( [(2 , 3), (2, 10), (1, 7), (1, 4), (7, 4), (4, 9), 

(5, 6), (6, 8), (3, 8), (3, 9), (8, 9)]) 



tionary: 




i" if i G nodes 
" if i £ nodes 



Fig. 2. The graphs X\ (left) and X 2 (right) in Example 1 



We test whether they are isomorphic: 

sage : Isomorphism2(X3,X4, 10, Iso=True) 

1 — > 2 

2 — > 1 

3 — > 7 

4 — > 4 

5 — > 5 

6 — > 6 

7 — > 3 

8 — > 8 

9 — > 9 

10 — > 10 
True 

And, indeed, it is obvious that this is an isomorphism between X\ and X<l. 
Example 2. Consider now the graphs X\ and X2 depicted in Fig. 3, 

sage: Xl=Graph( [(1 , 7), (1, 8), (1, 10), (2, 3), (3, 6), 
(4, 5), (5, 6), (6, 10), (7,9), (7, 10), (8, 9)]) 

sage: X2=Graph( [(1 , 7), (1, 9), (2, 3), (2, 5), (2, 10), 
(4, 5), (4, 6), (4, 10), (6,8), (7, 8), (7, 10)]) 

sage: Isomorphism (XI ,X2 , 10) 

False 

And, indeed, they are clearly non isomorphic. 

Test 2 Fig. 4 shows the time needed by Isomorphism2 to test whether two ran- 
dom graphs are isomorphic, as a function of the numbers of nodes in the graphs. 
The times are so small because two random graphs with the same number of 
nodes have probably different numbers of edges, a property that our program 
checks before proceeding with Luks' algorithm. Although in most cases the algo- 
rithm detected non-isomorphism by trivial reasons, the algorithm is also relatively 
efficient in non trivial cases. 



Fig. 3. The graphs X\ (left) and X 2 (right) in Example 2 




Fig. 4. Average time (in seconds) for different numbers of nodes when testing the 
isomorphism of random graphs 



Test 3 In this test we fix the degrees ofn—1 nodes in the graphs, and the degree 
of the last node is chosen at random. In this way we guarantee that the probability 
of two graphs being isomorphic is 1/3. In this case, our algorithm also runs in 
reasonable time: see Figs. 5 and 6. 



Fig. 5. Average time (in seconds) for different numbers of nodes when testing the 
isomorphism of semirandom graphs 



Fig. 6. The band under 2 seconds of the graphic in Fig. 5 



Test 4 Fig. 7 shows the times (red dots) needed by the algorithm to detect the 
isomorphism between pairs of isomorphic graphs with n nodes, and it compares 
this time with the functions (n/10) 4 ; (n/10) 3 ; (n/10) 2 log(n/10) ; and (n/10) 2 . 




Fig. 7. Comparison between the algorithm and the functions 
(n/10) 4 , (n/10) 3 , (n/10) 2 log(n/10), (n/10) 2 



Our last test deals with rooted phylogenetic networks. 

Test 5 Fig. 8 displays the relation running time-number of nodes for our algo- 
rithm when applied to random fully resolved rooted phylogenetic networks on the 
same sets of taxa, and Fig. 9 shows this relation for phylogenetic networks that 
are isomorphic as undirected graphs but need not be isomorphic as phylogenetic 
networks. 

4 Conclusions 

In this paper we have presented our implementation in SAGE of Luks' polynomial- 
time algorithm for testing the isomorphism of ternary graphs. This algorithm 
has interesting applications in phylogenetics, as it allows, for instance, to detect 




Fig. 8. Average time (in seconds) for different numbers of nodes when testing the 
isomorphism of fully resolved rooted phylogenetic networks 




Fig. 9. Average time (in seconds) for different numbers of nodes when testing the 
isomorphism of fully resolved rooted phylogenetic networks that are isomorphic as 
graphs 

whether two fully resolved rooted phylogenetic networks are isomorphic. There- 
fore, we have adapted and implemented Luks' algorithm for this type of graphs. 
Our adaptation has been, except for one point, a direct translation of Luks' algo- 
rithm. Fully resolved rooted phylogenetic networks have specific characteristics 
that could be used to improve the algorithm to make it still more efficient in this 
specific application. It is in our research agenda to develop such an adaptation, 
and we hope to present it elsewhere. 
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